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ABSTRACT 



Or 



Context. The Milagro experiment has announced the discovery of an excess flux of TeV cosmic rays from the general direction of the 
heliotail, also close to the Galactic anticenter. 

Aims. We investigate the hypothesis that the excess cosmic rays were produced in the SN explosion that gave birth to the Geminga 
j_( . pulsar. 

Methods. The assumptions underlying our proposed scenario are that the Geminga supernova occurred about 3.4 10 5 years ago (as 
indicated by the spin down timescale), that a burst of cosmic rays was injected with total energy ~ 10 49 erg (i.e., about 1% of a typical 
SN output), and that the Geminga pulsar was born with a positive radial velocity of 100-200 km s . 

Results. We find that our hypothesis is consistent with the available information. In a first variant (likely oversimplified), the cosmic 
rays have diffused according to the Bohm prescription (i.e., with a diffusion coefficient on the order of c x r L , with c the speed of light 
and r L the Larmor radius). An alternative scheme assumes that diffusion only occurred initially, and the final propagation to the Sun 
was a free streaming in a diverging magnetic field. 

Conclusions. If the observed cosmic ray excess does indeed arise from the Geminga SN explosion, the long-sought "smoking gun" 
connecting cosmic rays with supernovae would finally be at hand. It could be said that, while looking for the "smoking gun", we were 
CNl ' hit by the bullets themselves. 
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^ 1. Introduction 2. The heliospheric scenario 

r~j ■ In a recent paper dAbdo et al. 2008I) . the Milagro collaboration All known effects of the heliosphere on the cosmic rays (so- 
rN | reports the detection of an excess of cosmic rays from the general lar modulation, anomalous cosmi c rays) are detected at energie s 



direction of the heliotail, also close to the Galactic anticenter, of about 1 GeV per nucleon (e.g., |Ngobeni and Potgieter 2 008) 



The excess is diffuse but confined (in the following we assume a much lower than the several TeV observed by Milagro. On top 

solid angle of 0.3 steradians to account for both hot spots A and of this, one can give a direct counter argument based on the en- 

B), is composed of hadrons (photons and electrons are excluded ergy budget. In the following, we take data about the cosmic rays 

to a highly significant level) and has a harder spectrum than the from Longair (1981) and data about the Very Local InterStellar 

general cosmic ray background up to about 10 TeV. Medium (VLISM) from Axford and Suess (1994). 

The authors discuss some possible explanations and con- From the Milagro paper (their Fig. 4) one deduces that the 

elude that none of them is viable. In particular, because of the excess flux at 10 TeV -measured as a fraction of the background 

positional coincidence of the excess with the heliotail, they con- cosmlc ra y nux " amounts to 1 5 1 0" in region A and 6 1 0" in re- 

sider a local origin of the phenomenon, but discard it on the 8 lon B - These values refel " to the cores of the two regions, which 

ground that several-TeV particles could not be easily generated have a relatively small angular extent (for instance, the core of 

or confined by the heliosphere region A is only 0.02 steradians). To account for the lower level 

, ,.• v . • u ■ cut- u ■ ■ ■ excess that is visible around the cores, we add two times the 

In this Note, we revisit the hypothesis of a heliospheric origin . ,., , . „ _ ' . TT , 

of the hot spots and provide an additional quantitative argument core counts u ln a solld of , 03 Under these as- 

against it. Then we point out that the closest plausible extra solar s um P tlons the av f a g e fractional excess is 5 10~ 4 , and the excess 

... . ■ , . , A , ■ flux turns out to be 
source is the supernova that produced the Geminga pulsar, and 

show that there is a region in the parameter space where this . c ,„_4 „ ~ , _ 1A -6 , n 1fl -9 -2 -l „, 

. . 6 ,. , F F O ~ 5 10 x 0.3 x 6.7 10 ~ 1.0 10 erg cm s . (1) 
alternative hypothesis is valid. 

The acceleration region must be at least as large as one gyration 

Send offprint requests to: M. Salvati radius, which for a 10 TeV proton in a 1 //G magnetic field is 
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ri = 0.01 parsec. Given the angular extent of the excess, the dis- 
tance to the acceleration region must also be at least as much as 
ri , i.e., at least 20 times the distance to the heliopause. In the 
heliospheric hypothesis, the acceleration region is powered by 
the converging flows of the solar wake, and, at the inferred dis- 
tance, the convergence angle must be ~ 1/20 radians. Finally, the 
power dissipated in the acceleration region (for a local velocity v 
of 25 km s _1 and a local density n of 0. 1 cm" 3 ), is approximately 



P~d^n(^f 10 23 ergs" 1 . 

This corresponds to a maximum flux at Earth equal to 
P 



2.6 10" u ergcnT 2 s"\ 



Anr\ 



(2) 



(3) 



which is much too low with respect to the measured one. 



3. The supernova scenario, part one 

The basis of the supernova scenario is twofold. First, the Milagro 
excess flux comes from the right direction in the sky: hot spots 
A and B are about 50 degrees apart, and nicely encompass 
Geminga (e.g., Bignami and Caraveo 1996). The pulsar has a 
non-negligible proper motion (125 km s"" 1 at a distance of 155 
parsec) so that its position at birth was different from the present 
one, the more so if the birthplace was close to the Sun. For the 
cases of interest, the displacement is around 20-30 degrees to- 
wards the south of region A. Such angular distances do not seem 
implausible in view of the effect of the magnetic field on the 
arrival direction. Second, a 10 TeV proton diffusing in a 1 /jG 
magnetic field in the Bohm regime (i.e., D ~ c x r£) reaches 
an e-folding distance of 65 parsec in the time elapsed since the 



Geminga supernova explosion (t exp 
spin down age is adopted): 



3.4 10 5 yr, if the pulsar 



which suggests a velocity vector within 30 degrees from the 
plane of the sky. Based on this result, Pellizza et al. (2005) put a 
lower limit of 90 pc on the distance from the Sun at which the su- 
pernova explosion might have occurred. On the other hand, a ra- 
dial velocity as high as assumed here does not seem implausible 
with respect to the measured transverse velocity; moreover, an 
even higher value has been suggested (Geh rels and Chen 1 993) 
in an attempt to relate the Geminga supernova to the formation of 
the Local Bubble. In any case, the scenario proposed in the next 
section can accomodate a distance at the lower limit of Pellizza 
et al. (2005). 

The density distribution of particles diffusing with a constant 
coefficient in a 3-dimensional region is 



( ro< ) 



n(r, t) = N 



3 yTF ) 3/2 



(6) 



where N is the total number injected in a small volume at r — 
0, t — 0. If e is the energy of the particles, the net flux at radius r 
and time t is 



dn(r, t) , r 
<D = -eD — — = en(r, t)— . 

8r V '2t 



(7) 



If we set t = t exp , e — 10 TeV, and r — R and require agreement 
with the right hand side of Eq.(l), we deduce n from Eq.(7) and 
N from Eq.(6). Finally, we assume that the excess cosmic rays 
have the same spectrum as the background cosmic rays, and ob- 
tain the following estimate for the cosmic ray output E of the 
Geminga supernova 



n(R, t exp ; e) = 6.7 10~ 18 cm" 3 , N{e) = 1.7 10' 



-.45 



£ = 1.5 10 49 erg 



(8) 



(9) 



R = j4Dt^p = ^4cr L t exp = 65 pc. 



(4) 



The diffusion of cosmic rays in the Galaxy (e.g., 
IDe Marco et al. 20071 ) is described with two diffusion coeffi- 
cients, D|| and D ± , parallel and orthogonal to the magnetic field, 
respectively. The coefficients may be written as 



c x A, D ± ~ c x 



A 



(5) 



where A is the mean free path along the magnetic field. Usually 
A is taken much greater than ri, so that D\\ is much larger than 
D ± and the diffusion is strongly anisotropic. However, there are 
circumstances where the two become equal to each other and 
the diffusion is isotropic in the Bohm regime: this happens in the 
limiting case of a very chaotic magnetic field, with {6B/B) ~ 1 
over distances ~ ri, so that A ~ r^. It might be argued that the 
explosion of the Geminga supernova in the relatively recent past 
has stirred the local interstellar magnetic field up to the required 
level of chaos. We pursue this hypothesis as a zeroth order as- 
sumption, useful for simple calculations, and a more articulated 
scenario will be presented in the next section. 

The present distance to Geminga is estimated to be 
155*59 P c ' so we must assume a non negligible velocity of 
the pulsar in the positive radial direction, equal to at least 
160 km s _1 . On the one hand, such a radial velocity is discordant 
with the morphology of the Geminga trail (Cara veo et al. 2 003), 



This estimate is perfectly in line with the commonly required ef- 
ficiency (about 1 %) with which a supernova energy output must 
be channeled into cosmic rays if indeed supernovae are to main- 
tain the Galactic cosmic ray reservoir. 

The irregular distribution of the excess flux, and especially 
the presence of two disjoint hot spots, is perhaps a conse- 
quence of large-scale irregularities in the background medium 
and background magnetic field: the shape of the diffusing cloud 
must be much more complex than a perfect sphere. 

A final comment is in order about assuming a spectrum of the 
excess cosmic rays similar to the one of the background cosmic 
rays, while the Milagro data indicate a much flatter slope (1.5 
versus 2.6) and a cutoff above several TeV. The qualitative expla- 
nation that we propose has to do with the dependence of and 
D on the particle energy. If at about 10 TeV the e-folding point 
of the diffusing cloud profile has reached the Sun, at much higher 
(lower) energies the e-folding point is much beyond (before) the 
Sun position. The ratio of the excess flux at a generic energy to 
the one at the fiducial energy (10 TeV) can be expressed as a 
function of the ratio of the relevant diffusion coefficients D and 



/Diof 2 



(10) 



so that in both limits the excess flux is diminished with respect 
to the fiducial case. 
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4. The supernova scenario, part two 

Drury and Aharonian (2008) have raised two important objec- 
tions to the scheme presented above. First, a diffusing cloud of 
cosmic rays would produce a very wide signal in the sky, instead 
of the relatively narrow hot spots detected by Milagro. Second, 
while diffusion in the Bohm regime is thought to occur in pe- 
culiar regions, it cannot be the general process governing the 
propagation of cosmic rays across large distances. They suggest 
that the hot spots might be due to excess cosmic rays streaming 
almost freely from a magnetic nozzle along a diverging field; the 
source of excess cosmic rays should be relatively nearby (100 pc 
or less), located at the nozzle or behind it, and its energy content 
should be a fraction of a supernova output. 

We note that the cosmic rays cannot stream freely all the 
way from the source to the Sun, because the propagation time 
would be too short, less than a thousand years, and such a young 
supernova remnant could not remain unnoticed. Moreover, the 
scenario of Drury and Aharonian (2008) does not provide an ex- 
planation for the peculiar spectral shape of the Milagro signal. 
One must assume that propagation occurs by diffusion from the 
source to the magnetic nozzle, over a time long enough to al- 
low the dissipation of the SNR. This initial part of the propaga- 
tion process is very similar to what we discussed in the previous 
section. In particular, diffusion would again act as a "passband" 
filter in energy, producing a hard spectrum with a high energy 
cutoff, analogous to the observed one. In the new scheme we 
can relax the Bohm assumption A ~ r^, Dy ~ D ± , and can acco- 
modate a wider range of values for the distance of the supernova: 
indeed, at variance with the previous scheme, now the diffusing 
cloud of cosmic rays does not need to propagate from the su- 
pernova to the Sun, but only from the supernova to the magnetic 
nozzle. More precisely, diffusion is only needed from the rim of 
the supernova remnant to the magnetic nozzle, and across-held 
diffusion (which is the slowest process of all) is needed only 
from the rim of the SNR to the first "useful" magnetic line (see 
Fig. 1). 

We conclude by arguing that it is unlikely that the supernova 
responsible for the excess cosmic rays was not the same that 
produced Geminga. In a cone with vertex on the Sun, axis in 
the direction of Geminga, height 150, and base radius 50 pc, we 
expect less than 0.01 supernovae in 3.4 10 5 years for a Galactic 
rate of 0.01 yr _1 . 

In Fig. 1 we schematize the proposed geometry. The upper 
panel is the projection of the anticenter region on the plane that, 
being orthogonal to the Galactic plane, contains the present po- 
sitions of the Sun and Geminga (the latter at the "close" value 
of 120 pc). The lower panel is the projection of the same region 
on the Galactic plane itself. The directions of the Milagro hot 
spots are indicated, with the exception of hot spot B, which is 
not drawn in the lower panel since it is very wide in Galactic 
longitude and would cover all other elements. The two stars 
mark two possible positions of the supernova explosion, placed 
at distances of 90 and 65 pc from the present position of the 
Sun, respectively. The circles have radii of 10 pc, and indicate 
the volume occupied by fully developed supernova remnants. A 
possible magnetic nozzle is sketched with heavy lines: one sees 
that the cosmic rays, after leaving the remnant, need only diffuse 
across the field for very few parsecs before catching the right 
field line and propagating to the Sun. 

Region B is 50 degrees away from Region A and could have 
a different origin. However, this hypothesis would entail several 
additional ad hoc assumptions and would, in any case, conflict 
with the low probability of multiple nearby supernova explo- 




-150 -100 -50 

parsec 



Fig. 1. Projection of the anticenter region on the meridian plane 
at Galactic longitude 195° (upper panel) and on the Galactic 
plane (lower panel). See text for details. 

sions. It is still more economical to attribute both regions to the 
same source. Then the large angular separation should be as- 
cribed to large irregularities of the magnetic field in the Solar 
vicinity: some of the field lines diverging from the "primary" 
magnetic nozzle intersect the line of sight to region B, as shown 
in the upper panel, and could fuel a "secondary" magnetic nozzle 
in that position. Admittedly, this is another epicycle in the mo- 
del, but it would have to be invoked in any case, even if region B 
were fueled by a source independent of region A. 

5. Conclusions 

Our discussion does not prove that the Milagro anticenter hot 
spots are a transient relic of the explosion of the Geminga su- 
pernova. However, it does provide a consistent framework for an 
experimental result that would otherwise remain unexplained. 

If our proposed picture were indeed true, one would con- 
clude that after all supernovae do produce cosmic rays with an 
efficiency of about 1%. We could see this transient relic only 
because the supernova that gave birth to the Geminga pulsar 
exploded very nearby, and not very long ago. Because of such 
lucky circumstances, the long sought "smoking gun" connecting 
cosmic rays with supernovae would finally be at hand. Indeed, it 
could be said that while looking for the "smoking gun" we were 
hit by the bullets themselves. 
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